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Abstract. - In this paper we provide a new analysis of tlie system of partial differential equations 
describing the radial and vertical equilibria of the plasma in accretion disks. In particular, we show 
that the partial differential system can be separated once a definite, oscillatory (or hyperbolic) 
form for the radial dependence of the relevant physical quantities is assumed. The system is thus 
reduced to an ordinary differential system in the vertical dimensionless coordinate. The resulting 
equations can be integrated analytically in the limit of small magnetic pressure. We complete our 
analysis with a direct numerical integration of the more general case. The main result is that a 
ring- like density profile (i.e., radial oscillations in the mass density) can appear even in the limit 
of small magnetic pressure. 



Introduction. — The morphology of an accretion disk 
around an astrophysical compact object represents one of 
the most important open questions of stellar physics [Ij . In 
fact, while in absence of a significant magnetic field of the 
central object the disk configuration is properly described 
by the fiuid dynamics approach, the situation becomes 
puzzling when we deal with a strongly magnetized source, 
namely a pulsar, accreting from a less dense companion. 
As recently shown by B. Coppi, see [2]- [5], the plasma 
nature of the disk implies a significant coupling between 
the vertical and the radial equilibrium, as a consequence 
of the relevant Lorentz force acting inside the structure. 
The existence of such a coupling suggests a deep modi- 
fication of the original point of view at the base of our 
understanding of the stellar accretion phenomenon; for a 
sample of the basic literature in the field, see [HI- 0. In 
fact, the standard approach to the description of a thin 
disk relies on the idea that the vertical equilibrium can 
be averaged out when the viscoresistive MHD is applied 
to the plasma. Such a model seems to satisfactorily re- 
produce the coarse-grain phenomenology, but at the price 
of introdocing an anomalous resistivity of the disk plasma 
(unjustified by direct estimations), see Ref. |10[|llj . 
The analyses in Refs. [21 [3] demostrate that the details 



of the disk equilibria are relevant in establishing an os- 
cillatory local structure inside the disk. In particular, in 
Ref. |i3j it is shown that, for a disk having a sufficiently 
strong magnetic pressure, (i. e. a small enough ratio of 
the thermostatic pressure to the magnetic one), the mass 
density perturbations, due to the internal currents, are 
able to induce a ring-like profile. This ideal MHD result 
constitutes an opposite point of view with respect to the 
idea of a diffusive magnetic field within the disk, as dis- 
cussed in Ref. [T]. The striking interest in the details of 
the local disk morphology consists in the possibilty that 
jets of matter and radiation are emitted by virtue of the 
strong magnetic field and the axial symmetry, see for in- 
stance Ref. |12j . 

Here, we provide a novel analysis of the fundamental 
partial differential system derived in Refs. for the ra- 
dial and vertical equilibria in the disk plasma. Our study 
is based on a separable solution, able to reduce the coupled 
partial differential scheme to a simple ordinary differential 
system in the vertical dimensionless coordinate. The sep- 
aration is realized by a suitable trigonometric expression 
for the radial dependence; the remaining unknowns are 
four functions, associated to the vertical dependence of 
the magnetic flux function, the mass density corrections 
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and eventually to functions related to the different behav- 
ior of the thermostatic pressure term. 
Indeed it is just the presence of two different radial behav- 
iors of the pressure, the main new feature of the present 
analysis. We see that the mass density perturbations in- 
duced in the plasma are relevant even when the ratio of 
the thermostatic to the magnetic pressure is high (differ- 
ently from the analyis in [5]. This correction to the mass 
distribution has an oscillating character, so that we see 
the formation of the ring-like profile even in the parame- 
ter region where the relevance of the magnetic field is not 
crucial on the background level. This result suggests that 
the presence of rings in the disk structure morphology is 
a very general feature of the accretion disk structure for 
magnetized stars. 

The ordinary differential system we derive is analyti- 
cally integrated for small values of the ratio of the mag- 
netic pressure over the thermostatic one; the solution that 
is obtained is a good tool to fix the boundary conditions of 
the numerical analysis for the general case. In fact, the an- 
alytic solution remains valid for general values of the free 
parameters, as far as we restrict ourselves sufficiently close 
to the equatorial plane, where the boundary conditions for 
the numerical analysis can be given. The main implication 
of this link between the analytical and numerical analyses 
is that we get a direct relation between the ratio of the 
magnetic to thermostatic pressure to the one between the 
perturbation wavenumber and the fundamental wavenum- 
ber of the plasma structure. Then our solution cannot ex- 
plore the very extreme value of the parameters, where the 
magnetic pressure completely dominates the equilibrium 
configuration. The paper is organized as follows. In Sec. 
I, we describe the basic features of the disk. In Sec. II, we 
write the equations governing the radial and vertical equi- 
librium of the disk. In Sec. Ill we reduce the fundamental 
partial differential system to a system of ordinary differ- 
ential equations, and we solve it analytically in the limit 
of small magnetic pressure, discussing the appearance of 
an oscillatory structure. In Sec. IV we show the results of 
the numerical integration, and finally in Sec. V we draw 
our conclusions. 

Basic Features of the Disk. — The magnetic field, 
characterizing the central object, takes the form 



B = --dzipe,- + -6*0 + -dr^Cz , 
r r r 



(1) 



with ip — xp{r , z'^) and / — /(^ , z). 

The matter flux associated with the disk morphology is: 

ev = dzQcr + euj{r , z'^)re^ H — drQe^ , (2) 

r r 

where e denotes the matter density and Q{r , z) is an odd 
function of z, to deal with a non zero accretion rate, i.e. 



Md = -27rr / evrdz = ^TiQ{r , zq) = 2ttI > , (3) 



za{r) ^ r being the half-width of the thin disk. 

The similarity of the magnetic field and matter flux 
structure, is due to their common divergenceless nature. 
Since in the present analysis we are concentrating our at- 
tention on the formation of the ring profile within the disk, 
in what follows, we neglect the presence of the functions 
/ and Q, which are relevant for the characterization of 
the azimuthal equilibrium. In fact, as discussed in 
the origin of the oscillatory structure comes out by the 
coupling of the vertical and the radial equilibria, when 
the internal currents rising in the plasma are taken into 
account. 

We now develop a local model of the equilibrium, as 
settled down around a radius value r = tq, in order to in- 
vestigate analytically the effects induced on the disk pro- 
file by the electromagnetic reaction of the plasma. To this 
end we split the energy density and the pressure contri- 
butions as e = e(ro, z^) -I- e and p = p(ro, z^) -I- p, respec- 
tively. The same way, we express the magnetic surface 
function in the form il> — '0o(^o) + "01 (''o i ^ ~ ; -2^), with 
V"! <C f/'o- The quantities e, p and describe the change 
of the fundamental plasma functions due to the currents 
that emerge within the disk embedded into the external 
magnetic field of the central object. In general these cor- 
rections are small in amplitude but with a very short scale 
of variation. Thus, we are led to address the "drift order- 
ing" for the behavior of the gradient amplitude, i.e. the 
first order gradients of the perturbations are of zero-order, 
while the second order ones dominate. 

As ensured by the corotation theorem [13] , the angular 
frequency of the disk rotation has to be expressed via the 
magnetic flux function as uj{ip). As a consequence, in the 
present splitted scheme, we can take the decomposition 
Lu = iox + V-"! I where luk is the Keplerian term and 
lu'q = dujQdipo — const. This form for w holds locally, as 
far as (r — rg) remains a sufficiently small quantity, so that 
the dominant deviation from the Keplerian contribution is 
due to ipi. 

Accordingly to the drift ordering, the profile of the 
toroidal currents rising in the disk, takes the expression 
~ -(c/47rro) x + S^Vi). 

Vertical and Radial Equilibrium. We now fix 

the equations governing the vertical and the radial equilib- 
rium of the disk, by separating the basic fluid component 
from the presence of the electromagnetic reaction. Such 
a splitting of the MHD equations for the vertical force 
balance gives 



dzP 



1 



rr2 ^KbT 
^0 = 2~ 



(4) 



+ = 0, (5) 



where eo(ro) = e(ro, 0) and rrii is the ion mass. The be- 
havior of the function D{z'^) accounts for the pure ther- 
mostatic equilibrium holding in the disk when the vertical 
gravity (i. e. the Keplerian rotation) is large enough to 
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provide a confined thin configuration, while the tempera- 
ture T admits the representation 



2KbT. 



rrii 



P 



P + P 



2KB{T + f). (6) 



The radial equations underlying the equilibrium of the ro- 
tating layers of the disk, can be decomposed into the dom- 
inant character of the Keplerian angular velocity plus an 
equation describing the behavior of the deviation Suj: 



1 



47rr^ 



= dr 



P 



{dripi 



(7) 



We neglected, in the radial and vertical equilibria, the 
presence of the poloidal current associated with the az- 
imuthal component of the magnetic field. 

We define the dimensionless functions Y , D and P, in 
place oi ipi, e and p, i. e. 



(8) 



where po = 2i^sf eo/m^ and /3 = Sttpo/Bo^ = l/lSe^) = 
/cQ-ffQ/S. We introduced the fundamental wavenumber fco 
of the radial equilibrium, defined as fcp = 3a;|./u^, with 
v\ = Att€i^/B1q, recalling that Bzo — drgi/jo/rQ. It is then 
natural to deal with the dimensionless radial variable x = 
ko{r — tq) ^ while assuming that the fundamental length in 
the vertical direction is A = y^eji/o, leading to introduce 
u = z/A.l 

By this definitions, the vertical and radial equilibrium 
equations can be restated respectively as 



D^-D\Y- 



dl.Y + e^dl.Y + -d,P ^ 



+ {dl.Y + eX,Y)d^Y ^Q. 



(9) 



(10) 



Once D and D are assigned, the equations above provide 
a coupled system for P and Y , allowing to fix the disk 
configuration due to the toroidal currents. 

Reduction of the Fundamental System. The 

analysis of the partial differential system derived above has 
been performed in Ref. [2\ in the limit of small values of tz 
and an approximated solution was found as an expansion 
in such a parameter. Instead in Ref. [3J, the study has 
been extended to the case > 1, by requiring that the 
function Y satisfied the basic eigenstate equation 



dl.Y 



.di.Y = -DY. 



(11) 



Here we show that the two partial differential equations © 
and ([TU]) can be treated separating the radial and vertical 
dependence, thus reducing them to an ordinary differential 



system. In fact, we easily get such a reduction by the 
following positions 

F = sin(aa;) (12) 

P = L{u^)cos{ax) + M{u^)sin^{ax) (13) 
b = d{u^) cos{ax) . (14) 

and by imposing the vanishing of the coefficients of each 
type of trigonometrical terms. A simple calculation shows 
that the vertical equilibrium ^ yields the two equations 

dL , , 

(15) 



e,d = 



dM 
dv? 



d^F\ dF 



du^ 



dv? ' 



while the radial equation (fTOll gives 



D{u')F-a^F + ez'^-'^L = Q 



dF + aM + aF -a^F + e 



du^ 2 
d^F 



dv? 



= 0. 



(16) 

(17) 
(18) 



1 

These two pairs of equations form an ordinary differen- 
tial system in the variable u of four coupled second order 
equation in the four unknowns F{v?), L(v?), M(v?) and 
d(u^) respectively. The quantities and a (we recall that 
P = l/3el) play the role of free parameters of the problem. 
In particular measures the relevance of the electromag- 
netic interaction in the establishment of the equilibrium 
configuration of the disk plasma. The greater Cz is, the 
stronger the internal currents deform the background dis- 
tribution of matter and magnetic field. The parameter 
a fixes the amplitude of the radial wavenumber (with re- 
spect to the fundamental one fco) associated to the pertur- 
bations. The greater is a, the smaller is the wavelenght of 
the radial plasma structures. The same way, also can 
be regarded as the parameter which gives the scale of the 
vertcal confinament, according to the relation, introduced 
above, A = ^/e^Ho. Such a relation, together with the def- 
inition A = y^Ho/ko, allows to express the function D{u?) 
from D{z'^) introduced in (|4]) as D{u^) = cxp{— e^w^}. 

Analytical Solution for Small ez Values. Let us study 
the system of configuration equations in the limit of small 
values of the parameter ez , when we can use the expansion 
D{u^) 1 — e^u^. In this way, we are led to search a 
solution to the four equations above, in the form 

F = Aexp(-i^); L = IF; (19) 

C(u2)f2 . (20) 



kF: 



M 



In other words we assume that the function Y is confined 
around the equatorial plane and that the other functions 
can be expressed through F{u'^), in agreement to the struc- 
ture of the four equations. Substituting expressions ([T^ 
into the equations (fTTj) . we get the algebraic relations 



a 



L a ■ 



(21) 
(22) 
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Thus, for small values of e^, we are able to provide an an- 
alytic solution describing the detailed features of the disk 
plasma. We see that the wavenumber of the perturbations 
is not very different, in this limit, from k^, while the func- 
tion d is much greater than F and, as we shall see below, 
this is an important peculiar feature of this solution. 

We also remark that, since this solution relies on the 
expansion 6""^"=" ~ 1 — e^u , its range of validity is ac- 
tually broader than the ez ^ I region. In particular, the 
solution is still valid for ~ 1, provided that u <C l/i/eT 
(i.e., provided that we are close enough to the equatorial 
plane). The solution can also be continued in the > 1 
region by noting that in this case, according to Eq. (PT|) 
above, a would be purely immaginary, and the trigono- 
metric functions would become hyperbolic functions. We 
are then led to search a solution in the form: 

Y' = F'{u^)smh{a'x) (23) 

P' = L'iu^) cosh{a'x) + M'{u^) sinh^(a'x) (24) 

D' = d' (u^) cosh{a' x) , (25) 

where, as before: 

i^' = v4'exp(-2^) ; L' = l'F; (26) 
d' = k'F ; M' = C'{u^)F'^ . (27) 
Repeating the above procedure we find: 



(28) 



C'iu^)^-a'^ -e^u^ . (29) 



The Oscillatory Structure. Once the form of the solu- 
tion has been fixed, we can analyze the physical implica- 
tions for the disk structure. In particular, it is immediate 
to recognize that for < 1 the mass density distribution 
acquires the behavior 

eo P 

2 

= 1 - e^u^ + A{a€z)e^~ cos(ax) > . (30) 

We see that the perturbations to the mass density profile 
is an odd function of x and has an oscillating radial de- 
pendence. On the contrary, it is clear that the oscillating 
behaviour is not present in the > 1 regime, since in 
that case the density and pressure are expressed in terms 
of hyperbolic functions. We will then concentrate in the 
following in the analysis of the < 1 case. 
Since e^w^ is much smaller than unity, it is easy to realize 
that the positive character of the mass density is ensured 
by the request A < l/{aez)- The total pressure term in 
the disk plasma is given by 



P_ 

Pq 



ieiA 



-a cos(ax)-|- 



A (a^ - e,u^) e 



sin^ {ax) 



and the perturbation is an even function of x. The term 
in squared brackets can provide a negative contribution 
to the total pressure, but, in the limit of small enough 
<C 1, its weight is limited by the coefficient and the 
expression above, for values of u in the range of some units 
over the equatorial plane, rewrites as 



P_ 

Po 



^ ~ 1 -f 3eU 



1 / ^ 3 , 
— ;= cos(ax) -I- -Ae 
V3 2 



sm' 



'{ax) 



(31) 



(32) 

and the total plasma pressure is clearly positive when 
A < l/{ael). It can be seen that, since < 1, this is 
always ensured once the condition for the positiveness of 
the density {A < l/{aez) ) is fulfilled. However, in the 
general case, the positive character of p/pa in the point 
ax = TT, requires that 2ae^^exp[(e2 — l/2)w^] ^ 1, that 
constrains the possible values of A, with implications on 
the morphology of the ring profile. 

However, the very important feature we get, is that, 
differently from the analysis in Refs. [2l|3], here the mass 
density can have nodes even for small values of e^, i. e. for 
high f3 values of the plasma (see Fig. [Ij. This behavior 
is a consequence of dealing with a solution in which the 
quantity e is of order e^, instead of order like in 
We show the radial density and pressure profiles in Figs. 
[1] and [21 for different values of the vertical coordinate u, 
and for values of the parameters A — 100 and — 10^^ 
[the value of a is fixed by the relation (HI])]. It is clearly 
seen that there are "empty" regions, i.e. regions where 
f/eo *C 1, indicating the presence of a ring-like structure 
of the disk. This a general feature of the mass distribution 
described by Eq. (|30l) when AaCz — 1. The relevance 
of this result relies on the existence of a local ring profile 
in the disk plasma even if the magnetic pressure does not 
dominate the termostatic one. As a consequence, we can 
infer that the oscillatory structure of the disk is expect- 
edly a very diffuse phenomenon in accreting astrophysical 
sources. 

Comparison with previous works 

In this section we compare our results with those ob- 
tained in previous works, mainly focusing on Ref. |3j. The 
main assumption underlying the analysis presented in Ref. 
[3] is that Y satisfies the eigenvalue equation (fTTj). In 
our analysis, this equation is not satisfied given the po- 
sitions that we have made concerning the form the func- 
tions Y and P [see Eq. ([T2|l above[. In fact, substituting 
Y — F{u'^) sin {ax) in the above eigenvalue equation, we 
get: 

d'^F 

tz-r^-a^F + DF = Q. (33) 
dw^ 

On the other hand, we have that instead, in our analysis, 
the corresponding equation satisfied by F is Eq. ITT)) , 
namely: 

ez'^-a-F + DF^^L^,. (34) 

We recall that the function L{u'^) appearing in the addi- 
tional term is related to the cosine part of the pressure. It 
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A = 100, £^ = 0.01 




Figure 1: Top: Radial density profile of the disk for u = 0.5, 1, 2 
(black [solid], red [dashed] and blue [dotted[ curves, respec- 
tively) . Bottom: Normalized mass density as a function of the 
dimensionless radial and vertical coordinates x and u. In both 
panels, A = 100 and = 0.01. 



A = 100, €z =0.01 




2 4 6 8 10 12 14 



X 

A = 100, ez =0.01 




Figure 2: Top: Normalized pressure as a function of the di- 
mensionless radial coordinate x, for u — 0.5, 1, 2 (black ]solid], 
red ]dashed] and blue [dotted[ curves, respectively). Bottom: 
Normalized pressure as a function of the dimensionless radial 
and vertical coordinates x and u. In both panels, A = 100 and 
= 0.01. 



it is thus clear that, if L{v?) ^ 0, the two approaches will 
bear different results. So the basic difference between the 
two analyses can be traced back to the presence, in the 
pressure function, of a term proportional to cos {ax) . 

Numerical Analysis. — In the region of the param- 
eters where is not much smaller than unity, the system 
p3|) - P?)l has to be integrated numerically. To this pur- 
pose, we change the vertical variable to ^ = Observ- 
ing that = + ^ rewrite the configuration 
system in the form: 



that, when ^ is small enough, the analytic solution found 
in the previous section is an acceptable solution to the 
differential system even if ez is of order unity or larger. 
In fact, the parameter of the expansion that leads to the 
above analytical solution is e^w^ = 2ez£,, so that for every 
value of Sz it is always possible to find a (however small) 
starting value for the numerical integration such that 
2ez^o < 1- 

Then, we give the initial conditions so that they match 
the analytic solution in 



^ + 2ezd = (35) 



dM 
~d^ 



2C 



d^ 



e-^^'^F-a^F + ez 



dF 

^d'^F dF\ 



dF 
~d^ 

a 
" 2 



SetdF + aM + aF 



^<PF_ 

He 



dF 



. (36) 
L = (37) 
= 0. (38) 



Fq = F(Co) = Ae-«° ; L{^o) - ^aF^ ; 
rf(?o) = ^Fo- Af(eo) = {a" - 2ez^o)Fl, 



dF 



(39) 



This system of four ordinary differential equations in the 
four unknown functions F(^), L(^), (^) and d{£,) can be 
integrated numerically by fairly standard methods. How- 
ever, in order to perform the numerical integration, initial 
conditions should be given. These can be found by noting 



The fact that we are matching the analytical solution 
in ^0 means that we have to enforce the condition a — 

-\/3(l — ez) like we did before, in order to ensure the 

consistency of the differential system. Hence, initial con- 
ditions are completely fixed once A and (or equivalently 
a) are fixed. The results of our numerical integration are 
shown in Figs. [3]|4l 
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A = 1 , £^ =0.9 




Figure 3: Behaviour of the functions F{u), L(u), M(u) and 
d{u) for D{u) = e''="^ A = 1 and = 0.9 (Red [sohd] Unes). 
Also shown is the correspoding analytical solution for D{u) — 
1 — EzM^ (Black [dotted] lines). 

Concluding Remarks. — In this paper we presented 
an exact separable solution of the radial and vertical equi- 
librium equations at the basis of the oscillating morphol- 
ogy emerging in the plasma configuration of an accretion 
disk. We reduced the original system of two partial dif- 
ferential equations into a set of four independent ordinary 
differential equations in the vertical coordinate, one for 
the flux surface function, one for the mass density per- 
turbations and two related to different components of the 
pressure. This reduced system admits an analytical solu- 
tion only in the limit of high values of the /3 parameter of 
the plasma, i. e. only when the thermostatic pressure of 
the disk is much larger than the magnetic pressure con- 
tribution. We also performed a numerical analysis of the 
system for small values of /3. 

The main results of our analysis can be summarized in 
the following three points. 

i) -We have derived a solution showing how the radial 
gradient of the thermostatic pressure is relevant in estab- 
lishing the equilibrium, even in the linear limit and for 
small values of e^. This feature makes our solution in- 
trinsically different from the analyses developed in Refs. 
[2l[3], except for the discussion of the extreme non- linear 
regime in Section IX of Ref. [3]. In fact, in such a limit, 
our approach is reconcilied with that one, because the lat- 
ter accounts on an equivalent level the role of the radial 
gradient of the pressure. 

ii) -We have obtained the fundamental feature that the 
ring-like structure emerges as a strong structural feature 
of the plasma disk confined in magnetic field, since the 
radial oscillation of the mass density takes place even in 
the linear weakly magnetized limit A <C 1 and <C 1. 
Indeed a suitable choice of these free parameters of the 
model is always possible in order to arrange for nodes in 
the mass density profile. This morphology is relevant be- 
cause suggests that the ring profile can be expected to be 
a general character of the magnetized accreting compact 
objects we observe in the Universe, and it is due to the 
direct link existing between the radial pressure gradient 
and the mass density perturbations. This output of our 



A= 1 ,e^ = 1.5 
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Figure 4: The same as Fig. |4l but for — 1.5 



separation algorithm indicates that more general configu- 
ration scenarios can be contained in the radial and vertical 
equilibrium equation with respect to the one investigated 
in though they could exist in a non-separable regime. 

iii)-We have found that the disk profile undergoes a 
transition when the parameter becomes greater than 
one, going from the oscillating structure to an hyperbolic 
behavior. This fact suggests the possible co-existence of 
two different disk components in the same global profile. 
In fact in the present local model we addressed as a 
constant because it refers to a generic value of the radial 
coordinate Tq, but throught the disck it is clearly a func- 
tion of the radial coordinate itself, taking values above and 
below unity in different radial regions. 
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